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Abstract: We present a study of loffe times in deep inelastic electron-proton scattering. 
We deduce 'experimental' loffe-time distributions from the small-x HERA data as described 
by a particular colour-dipole-model fit. We show distributions for three representative j*p 
cm. energies W and various values of the photon virtuality Q^. These distributions are 
rather broad for transversely and very narrow for longitudinally polarised virtual photons. 
The loffe times for W = 150 GeV, for example, range from around 10^ fm for = 1 GeV^ 
to around 10 fm for = 100 GeV^. Based on our results we discuss consequences for the 
limitations of applicability of the dipole picture. 
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1. Introduction 

In this article we shall present a quantitative study of the space-time structure of deep 
inelastic electron- and positron-proton scattering (DIS) 

e^+p^e^ + X (1.1) 

as measured extensively in HERA experiments |l], ^, ^, 0, |5|. We shall be concerned with 
the kinematic region where only exchange of a virtual photon between the leptons and the 



hadrons matters. Thus, as usual, with the reaction (IJ.) we study in fact the absorption 
of a virtual photon on the proton 

j*+p^X. (1.2) 

The total cross section for ( |1.2| ) is, apart from kinematic factors, given by the imaginary 
part of the amplitude for forward Compton scattering of a virtual photon on the proton, 

7*+P^7*+P- (1-3) 

The study of the space-time structure of this process has a long history going back to 
the discussions of the vector-meson-dominance (VMD) model in the 1960s and in particular 
to articles by Gribov et al. Q and loffe Q . In Q the question was posed where and when 
the initial virtual photon is absorbed, that is, where it fluctuates into 'hadronic stuff' and 
where the 'hadronic stuff' fluctuates back to give the final virtual photon. In this article we 
shall give quantitative answers to this question, as far as the region of small x is concerned, 
using the HERA data. In fact we shall use a fit to these data in a particular colour dipole 
model. 

Let us first make some more historical remarks. It is an old idea that the scattering of 
a highly-energetic photon on a hadron may essentially be considered as a strong interaction 
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process where the photon acts in some way as a hadron. Vector mesons as intermediary 
particles in the couphng of photons to nucleons were introduced in ^, |l^] in order to 
understand properties of the electromagnetic form factors. An important role was assigned 
to vector mesons in strong interactions in |11]. Then in [12| the vector-meson-dominance 
model was proposed. There it is assumed that whenever a photon couples to hadrons it 
first converts to the vector mesons p, uj, (j) (proposed / known at the time) with universal 
coupling constants. These vector mesons have then normal hadronic reactions. The VMD 
model was rather successful in describing many results from photon-hadron reactions. For 
reviews see, for instance, [|l3|, But the results from deep inelastic scattering presented 
problems. Indeed, in [15| a simple form of this VMD picture was applied to DIS and 
predictions were made which, however, were subsequently disproven by experiments. This 



led to the formulation of the concept of generalized vector meson dominance, see 14, 17| 
and references therein. The key idea in these descriptions is that the photon fluctuates 
into a series of vector mesons which subsequently scatter off the proton. Quite obviously, 
the lifetime of the fluctuation of the photon into a vector meson needs to be sufficiently 
long for that picture to be consistent. 

The colour dipole model ^] builds on similar ideas but is motivated to a large 

extent by perturbative QCD. In this picture the reaction (L2) is viewed as a two-step 
process. In the first step the photon splits into a quark-antiquark pair which represents 
the colour dipole. Subsequently that pair scatters on the proton, this second step being 
a purely hadronic reaction. For the applicability of the dipole model it is crucial that the 
lifetime of the fiuctuation of the photon into the colour dipole is much larger than the 
typical timescale of the dipole-proton interaction. In the context of the dipole picture the 
lifetime of the dipole fluctuation is usually referred to as the loffe time. We may note in 
parentheses, however, that loffe in his original paper Q was concerned with a somewhat 
different time as we shall comment on below. It is the aim of the present paper to study 
the distribution of loffe times and their dependence on the kinematical parameters of the 
photon-proton scattering process using the HERA data on DIS. 

The splitting of a photon of high enough virtuality into a quark-antiquark pair 
can be calculated in perturbation theory. The subsequent scattering of the colour dipole 
off the proton, on the other hand, is a genuinely nonperturbative process. Therefore, this 
second step of the 7*p scattering process needs to be described by suitable models. A 



variety of such models has been constructed, see |21, 23| for some prominent examples 



and 1 24, 25 1 for overviews. These models are very successful in describing the structure 
functions measured at HERA. 

But despite the impressive phenomenological success of dipole models there are some 
caveats. In p6|, 27 1 the foundations of the dipole picture were examined in detail. It 



was shown that a number of assumptions and approximations is required to arrive at the 
dipole picture. These results naturally raise the question about the range of validity of 
these approximations and assumptions. In ^] it was found that already the general 
formulae of the standard dipole approach allow one to derive stringent bounds on various 
ratios of structure functions. These bounds were then used to determine the kinematic 
region where the dipole picture is compatible with the data. An important point in this 
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connection concerns the variables on which the dipole-proton cross section depends. The 
natural variables were found to be r, the transverse size of the dipole, and W, the cm. 
energy of the dipole-proton scattering. Using this functional dependence for the dipole- 
proton cross section we obtained the following result: For 7*p c. m. energies W in the 
range of 60 to 240 GeV the standard dipole picture fails to be compatible with the HERA 
data for photon virtualities larger than about 100 to 200 GeV^ ||2^. Clearly, these 
large virtualities correspond to relatively short loffe times. In the present paper we want 
to study this relation between loffe times and the kinematical parameters in more detail 
by calculating the distribution of loffe times and the corresponding contributions to the 
structure functions for given values of and W. 

Let us, therefore, consider the reaction ( |1.2| ) at high c. m. energy in the proton rest 
system. The virtual photon 7* of 4- momentum q = {q^,q) fluctuates into a quark and 
antiquark of 4-momenta k = {k^,k) and k' = {k'^,k') respectively, where 3- momentum 
is conserved, q = k + k'. Energy, however, is not conserved. Instead, there is an energy 
mismatch 

AE = k^ + k'^ - g° (1.4) 

between the quark-antiquark pair and the virtual photon. According to the uncertainty 
relation such a fluctuation cannot live longer than the time 

(1.5) 



AE fco + fe'O - ■ 

In the following we shall call r the loffe time for the initial 7*. Discussions of loffe time 
distributions in the context of the operator product expansion can be found in |3C]. In 



1 31] the loffe-time structure of the gluon distribution function in the double logarithmic 
approximation was considered. 

In applications of the dipole model one frequently finds simple estimates for loffe times, 
and if they turn out to be of the order of several femtometers or larger this is taken as 
justification for using the dipole model. But the actual situation is more involved. We shall 



find that even for a fixed kinematical point for reaction (1.2) one has a whole distribution 



of loffe times. In fact, the 7*p total absorption cross section is most conveniently obtained 
from the imaginary part of the forward scattering amplitude for the reaction ( |1.3D , j*p — > 
7*p. As a consequence, we shall have to deal with two loffe times, one for the initial state 
7* and one for the final state 7*. Both times have distributions which also depend on the 
polarisation, transverse or longitudinal, of the 7*. We shall calculate such distributions in 
the following from the HERA data using a phenomenological dipole model which describes 
the data quite well. We shall choose the model of Golec-Biernat and Wiisthoff (GBW) |^ 
for this purpose. 

In section |2| we review the relevant formulae of the dipole picture and define the AE 
and loffe time r distributions in this framework. In section ^ we recall the basics of the 
GBW model for the dipole cross section. We then present numerical results for the AE 
and the corresponding loffe-time distributions. Our conclusions are drawn in section ^. In 
two appendices we present the details of our calculations. 
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2. loffe time and distributions in the dipole picture 



We use the standard formulae for the kinematics and for the definitions of structure func- 
tions of deep inelastic electron- and positron-proton scattering (|l.lD , see for instance p2|| . 
For momentum transfers squared < 1000 GeV^ it is sufficient to consider only the ex- 
change of a virtual photon between the lepton and the proton in ( |1.1| ). Thus, the reaction 
which we shall study in the following is the absorption of a virtual photon 7* on the proton, 

l*{q)+p{p)^X{p'), (2.1) 

where we indicate the 4-momenta in brackets. The cm. energy for this reaction is denoted 
by W , the virtuality of the 7* by Q^. For these and for Bjorken's scaling variable x we 
have 



W = {p + qf , 

^' = -^'' (2.2) 

g2 q2 



2p-q W"^ + -mj 



The proton in ( |2.1[ ) is supposed to be unpolarised, but the virtual photon can have trans- 
verse or longitudinal polarisation. The corresponding total cross sections are aT{W,Q'^) 
and a^iW^Q"^), respectively. The structure function F2 is 

F2{W,Q') = ^^[aT{W,Q')+aUW,Q')]{l-x) + o(^^^ . (2.3) 
For small Bjorken-x, x <C 1, this simplifies to 

F2{W, Q2) = [aT{W, Q2) + aiiW, Q^)] . (2.4) 

In the following we shall use this simpler relation since we shall only consider data for x <C 1. 
In order to obtain the standard dipole model for the cross sections ctt.l we can relate them 
first to the imaginary part of the virtual Compton forward scattering amplitude, 

7*iq)+pip) ^i*iq)+pip)- (2.5) 

The latter is then represented as the initial 7* splitting into a qq pair, this pair scattering 
on the proton, and the qq subsequently fusing into the final state 7*, see figure |^. Note 
that this figure is to be read from right to left in order to be in complete analogy with 
the occurrence of the various factors in the amplitudes in the equations below. With the 
assumptions spelled out in detail in section 6 of the diagram of figure |l] gives 

aTiW,Q^) = Y^ I dV rdaj;(^S±(a,r,Q))*<7W(r,H^)V'S,^(a,r,Q), (2.6) 

q A,A' 

^l(H^,Q') = E / rdaj;(^g,^(a,r,Q))*a(^)(r,W)V'S''(a,r,Q). (2.7) 

q A,A' 



- 4 - 



I 



Figure 1: Basic diagram for the description of the cross sections <jt,l of 7*75 scattering in the 
standard dipole approach 

Here a is the longitudinal momentum fraction of the 7* carried by the quark, r is the 
vector in transverse position space from the antiquark to the quark, r = |r|, and A and 
A' are the helicities of q and q, respectively. The total cross section for the scattering 
of the qq pair on the proton is denoted by a'^'^\ the wave functions for transversely and 
longitudinally polarised virtual photons by V'a'a'^ ^aa'^' ^respectively. For the explicit 
form of these wave functions see appendix A sum over all contributing quark flavours 
g is to be performed. 

However, the standard representation of the photon wave functions in terms of lon- 
gitudinal momentum fraction a and transverse position r is not suitable for a discussion 
of loffe-time distributions. To study them we have to go to longitudinal and transverse 
momentum space. There we can directly read off the energy mismatch between the 7* and 
the qq pair, both for the initial and the final states. Let us, therefore, consider first the 
splitting of the initial 7* into a qq pair in momentum space. We have 3-momentum but not 
energy conservation at this splitting. Taking this into account we define the 4-momentum 
of the quark, k, and that of the antiquark, k' , by 

k = aq + kx , k^ = \J k^ + , 

k' = {l-a)q-kT, k'° = -^fe'^ + ml . (2.8) 

Here, in essence, a £ [0,1]. The precise a-range is discussed in appendix ^ Then the 
photon wave functions in momentum space are easily calculated. For a derivation see for 
example p7[ ] . We give the results in appendix ^ both in the - in leading order in as and 
Oem ~ exact and in the approximate form which is usually used in the dipole model fits. 
The approximation involves in particular neglecting terms such as (fcy + mg)/(a^q^) and 
{k"^ + rrig) / {{1 — a)^q^) with respect to 1. However, for some given |g| those terms become 
non- negligible if fcy = I^tI is large or if a is close to or 1. If relevant contributions 
to some observable depend on the photon wave functions in this kinematical region the 
above approximation could become invalid. This is potentially important when considering 
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distributions in the loffe times, in particular for short loffe times. We shall calculate the 
distributions with and without the above approximation in order to quantify this effect. 

The photon wave functions in transverse position space are related by a Fourier trans- 
formation to their momentum space representations 

^S,±'^(a, r,Q) = l^ e''=-^V^S,±'^(a, kr, Q) . (2.9) 



The next step is to insert the representation ( |2.9D for the wave functions for both the initial 
and final state photon into ( p.6[ ) and (2.7). This gives 

i2Z j2; 



(2.10) 

aLiW,Q^)= ^ / d« '^^y (^2$ (ffA'''(«> fcT, Q)) * (kr - kr, W^)V^g,^a, kr, Q) . 

(2.11) 

Here we denote the Fourier transform of the dipole-proton cross section by 

~a^'i'^ [kr, W)= [ dV e**'^^ a^"?) (r, W) . (2.12) 



The four-momenta of the quark and antiquark in the initial state dipole are given in 
( |2.8| ). For the quark and antiquark in the final state dipole we denote the 4-momenta by k 
and k' , respectively. They are obtained from ( p. 81 ) with the replacements k ^ k, k' ^ k' , 
and fcr kx- Note that a stays the same. 

We recall now the definition of the energy mismatches AE for the initial 7* splitting 
to qq and AE for the final qq fusing to 7*: 

AE = k^ + k'^-q\ (2.13) 
AE = F + k'^-q°. (2.14) 

The corresponding loffe times are 

(2.15) 



AE 

f = ~- (2.16) 

We have AE > and AE > which implies also r > and f > 0. 

From ([2T0D and (|T1|) we see that the cross sections ax and ai involve the super- 
positions of amplitudes of various AE in the initial and various AE in the final state. 
Therefore, we define the joint distributions in 77 = AE and fj = AE by 

d''aT,L{W,Q^V,v) fd^kr f d^kr 



X a^'i^kT - kr, W) 6{7] - AE) V^S,^'^(a, kr, Q) . 



(2.17) 
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Note that these distributions are real- valued (see ( |A.8D ) but for rj fj they need not 
be positive. For specific dipole models we can now calculate a^*^^ and evaluate the two- 
dimensional distributions ( p. 17 ). However, one-dimensional distributions are more easily 



visualised. Thus, in the following we shall study how the sum of AE and AE 

AE+ = AE + AE (2.18) 

is distributed. We, therefore, define 

daT,L{W,Q^AE+) r^_ d^aT,L{W,Q^V,v) 

Me; = i ^^io ' d^v ^(^^^-v-v) (2.19) 

and the lofFe time corresponding to AE^ by 

1 TT 

(2.20) 



A^; I T + f 



The quantity 2r+ is called the 'harmonic mean' of r and r [33|. From ( p. 20 ) we get then 



r > r+ and f > r+. That is, for given r+ only individual loffe times r and f which are 
larger or equal to r+ contribute. At this point we may also note that loffe in his original 
paper considered the time between the initial 7* fluctuating into 'hadronic stuff' and 
this 'stuff' fluctuating back to the final 7*. In our calculation this corresponds to the time 
T + f. For given r+ we have t + f > 4:T^. 

In the next section we shall present numerical results for normalised distributions in 
lnA£'+ respectively lnr+ = — InAE'^. First we shall study these distributions for F2, 

AE+ dF2 _ AE+ f dar da^ \ 
F2 dAE+ aT + (TL\dAE+ dAE+J ' ^' ' 

It is also of interest to study the distributions in the loffe time r+ and in AE^ for the 
contributions of individual quark flavours to gt and cr^. For this we define and 



d'^afj^l dr]drj as in ( |2.10| ), ( 2.11 ), and ( |2.17] ), respectively, but omitting the sum over the 

quark flavours q on the r.h.s. Then we define daf^/dAE^ in analogy to ( |2.19| ). The 
corresponding normalised distributions in In AEj^ are given by 

AE+ daf]^ 



M) dAE, 



(2.22) 



3. Results 



In this section we present results for the distributions in AEj^ and in loffe times t+ in 
DIS. To evaluate ( |2.21| ) and ( |2.22D we have to choose a specific dipole-model fit to the 
data since we need the values for the dipole-proton cross sections a'^'^^; see (^.61) , (^3), and 



( 2.12 ). In the following we choose to work with the model constructed by Golec-Biernat and 



Wiisthoff ||2l|. This GBW model describes the F2 data from HERA quite well. Whether 



it also describes the longitudinal structure function Fl is not yet clear since the first Fl 
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measurements from HERA [34| are not precise enough to draw firm conlusions. The dipole 
cross section of this model is given for quark flavour q by 



where 



1 — exp 



4r^(i) 



(3.1) 



^o(^) = 7^(-) , 5^ = ^(1 + ^^ I- (3-2) 



We choose their parameter set which includes the charm quark, Qq = 1 GeV, o"o = 
29.12 mb, A = 0.277, xq = 0.41 • 10"^, nig = 0.14 GeV for q = u,d, s, and = 1.5 GeV. 
The 6-quark contributions are neglected in this fit. 

The GBW dipole cross section <5^Qg^ depends on r and Bjorken-x, and therefore not 



only on r and W but also on Q , see (|2.2D . As was discussed in ||2^, ^ the natural and 
- in our opinion - correct energy variable for the dipole cross section a^'^^ is W, and a^"^^ 
should actually be independent of Q^. A dependence on x requires additional assumptions 
which appear difficult to justify. We will discuss the problem of choosing the correct energy 
variable in more detail elsewhere. For the present considerations, however, this issue how 
the a^'^^ depend on the kinematic variables is not of prime relevance since we only need a 
dipole model fit to the data at given values of these variables. 

Now it is relatively straightforward to calculate the Fourier transform of <5"qbyv' 
sert it as well as the expressions for the 7* wave functions from appendix^ in ( 2.10 ), ( 2.11| ) 
and ( 2.17 ), and to obtain the numerical results for the distributions in InAE^ = — lnr+. 
The calculational details are presented in appendix The results for the normalised F2 
distribution ( 2.21 ) are shown for 7*p cm. energies W = 70, 150 and 220 GeV and a 
number of values in figure 0. In figure we present the normalised distributions in 
In AE+ = — In T+ for a^! and for the individual quark fiavours separately, here for the 
energy W = 150 GeV. 

In the calculations for figures ^ and ^ the full expressions ( |A.1| ) and (A. 2) are used 
for the photon wave functions rather than their high energy and small-A:^ approximation. 
Also, we integrate over the full range in a which is slightly bigger than [0, 1], see ( |B.ll ). 
We have studied the effect of the high energy approximation ( A.10| )- (|A.13| ) for the photon 
wave functions as well as of restricting the a integration range to [0,1]. The loffe-time 
distributions in figures ^ and |3| are essentially unaltered when performing the a integration 
only over [0, 1]. The a range specified by (|B.ll ) is only slightly larger than [0, 1] and no 
sufficiently strong enhancement of the integrand compensates for this fact in the kinemat- 
ical ranges considered here. As expected, effects in the loffe-time distributions from the 
high energy and small-A:^ approximation for the photon wave functions become larger for 
increasing AE-^-. These deviations would alter the curves of figures || and ^ to an amount 
that is clearly visible only in the high and small W cases. However, for all distribu- 
tions they are still small enough to be safely omitted from our further discussions here. We 
stress that this result is not obvious, since the interplay of the high energy and small-fcy ap- 
proximation and the dipole energy mismatch is non-trivial in the longitudinal-momentum 
endpoint regions. 
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Figure 2: The normalised F2 distribution (2.21) in the joint dipole-energy mismatch AE+ for the 



GBW model with charm. Supplementary to the AE+ values the corresponding values for the loffe 
times T+ are denoted for the abscissa. 

Let us now discuss the main features of the curves in figures |2| and ^. We first note 
that due to the normalisation of the curves in figure ^ one cannot deduce from them the 
relative importance of hght quarks and the c quark in F2. This information is obtained 
from figure where the secondary peaks in the individual curves are due to the c-quark 
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Figure 3: The normalised cr^'' and ct^'' distributions ( 2.22 ) in tlie joint dipole-energy mismatch 



AEjf. respectively loffe time for the GBW model with charm. 

contributions. The peaks seen in figures ^ and ^ are not infinitely sharp. We analysed 
their structure using an enlarged AEj^ scale and find them to be smooth finite maxima in 
all cases. 

The curves in figure |3| clearly exhibit the respective thresholds due to the kinematical 
lower bounds for the various quark flavours q, 



Ai?+,^in = 2 Jg2 + 4m2 - = 2 Jq^ + Am^„ - Vq^-Q^) , (3-3) 



see ( |A.20| ) and (2.18). The quantities A£^_|_^min are shown in the left graph of figure]^ for 



light quarks as functions of for W = 70, 150 and 220 GeV. In the right graph of figure ^ 
we show Ai?-|-^inin as a function of for W = 150 GeV for the light quark mass and for the 
charm mass. From figure ^ we see that the AE^ distributions for transverse 7* polarisation 
are significantly broader than those for longitudinal 7* polarisation. For both polarisation 
types the distributions become narrower for increasing values of Q^. The distributions for 
longitudinal polarisation have negative parts which is alright, see the discussion after ( 2.17] ). 
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Figure 4: Minimal energy mismatch AiJ+^min and corresponding maximal lofFe time r+.max of the 
7*-dipole transition as a function of the photon virtuality Q^. The curves are for different values 
of the energy W in the case of light quarks (left graph) and for different quark masses in the case 
= 150 GeV (right graph). 



The normalised distributions in T4. respectively AE^ shown in figure |2| are obtained 
from the DIS small-x data from HERA, using a specific dipole-model fit. In order to test 
the model dependence we have also investigated another fit of the GBW type and find 
compatible results. 

From figure § we can clearly see that the loffe times decrease with increasing at fixed 
W. At fixed the lofFe times increase with increasing W. The distributions are always 
rather broad. At = 70 GeV and = 100 GeV^, for example, the upper cutoff for the 
loffe time r+ is 5 fm and, indeed, the peaks of the distribution are there. But we have also 
significant contributions for t+ down to 1 fm. As stressed already in the introduction, for 
the dipole picture to be a valid description of DIS the loffe times have to be appreciably 
larger than the typical interaction time Thad of the dipole with the proton. As an estimate 
of the latter we take the electromagnetic radius of the proton Thad ~ ~ 1 fm. A necessary 
condition for both dipole lifetimes r and f (see ( 2.15D , ( |2.16| )) to be much larger than r^ad 
is thus A£^_|_ <^ 1/Thad or, equivalently, 

T+ > Thad ~ 1 fm . (3.4) 



If this condition is violated for a relevant portion of the r+ distribution of F2 the applica- 
bility of the dipole model is questionable for the kinematical point under consideration. 

From figure |2| we see that for I^ = 70 GeV the typical loffe times r+ are reasonably 
above 1 fm for all < 100 GeV^. There, the necessary separation condition ( p.4| ) is 
satisfied. However, for = 1000 GeV^ we see that the separation condition ( |3.4| ) is 
clearly violated, since the complete r+ distribution lies at t+ < 1 fm. For W = 150 GeV 
the distributions are shifted to larger values of r+. But still for = 1000 GeV^ the t+ 
distribution only starts at t+ ~ 2.5 fm and extends well below t+ = 1 fm. For W = 220 GeV 
and = 1000 GeV^ the situation is similar with the maximal r+ ~ 5 fm. Therefore, we 
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must conclude from our study of loffe-time distributions in DIS that the basic separation 
condition of time scales (3.4) is violated in the HERA energy range = 70 to 220 GeV, 
for values of hundred to several hundred GeV^ . 



4. Conclusions 



We have calculated loffe-time distributions for DIS in the HERA energy range in the 
framework of the dipole model. As a convenient fit to the data we used the Golec-Biernat- 
Wiisthoff parametrisation of the dipole-proton cross section. We have obtained loffe-time 
distributions for 7*p c. m. energies W = 70, 150 and 220 GeV and values ranging from 
0.1 to 1000 GeV^. A basic requirement for the dipole model to make sense is that the 
loffe times at the kinematic point considered are much larger than the hadronic timescale 
7"had ~ 1 fin- We find that typical loffe times are large with respect to this hadronic scale 
for Q2 < 100 GeV^ such that no inconsistencies arise there. However, at photon virtuality 
Q2 = 1000 GeV^ typical loffe times are of similar order as or smaller than the hadronic 
scale, which violates the standard assumption for the validity of the dipole picture. Thus, 
we find that in the HERA energy range the dipole picture starts to lack physical justification 
for values in the hundred to several hundred GeV^ range. 

Note that we come to this conclusion here using a GBW fit to the HERA data where 
the dipole-proton cross section is assumed to depend on r and Bjorken x. In [^] we have 
investigated the limits of applicability of the dipole model to the HERA data using the - in 
our opinion correct - dependence of the dipole-proton cross section on r and W. From figure 
9 of 1 29 1 we see that, nevertheless, we found a very similar range for the applicability of 
the dipole model. The upper limits for as obtained from this figure range from about 
120 GeV^ at ly = 70 GeV to about 180 GeV^ atW = 220 GeV. These are strict bounds in 
the sense that the data cannot be fitted for a larger range with a non-negative dipole- 
proton cross section depending on r and W. In our present work we find that even if a fit 
to the data with a dipole-proton cross section depending on r and Bjorken x is possible 
also for higher values of its physical meaning may become questionable there since the 
loffe times become too short. For very small values, < 2 GeV^ say, the loffe times 
shown in figure ^ are very large. Nevertheless, this does not immediately imply that the 
standard dipole model is without problems there. As discussed in |27, |2^, |2^ the lowest 
order expressions for the photon wave functions are expected to become unreliable in this 
kinematic region. 

In summary, we have calculated 'experimental' loffe-time distributions from the small-x 
HERA data as described by a GBW dipole-model fit. We have studied their dependence 
on the energy W and on the photon virtuality Q^. The loffe-time distributions of the cross 
sections are found to be rather broad for transversely and very narrow for longitudinally 
polarised virtual photons. Accordingly, the loffe-time distributions of F2 are always rather 
broad. 
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A. Photon wave functions and energy mismatch AE 

In this appendix we collect the formulae of the dipole model which are relevant for our 
calculations. The assumptions needed to arrive at the standard dipole picture of DIS are 



spelled out in detail in section 6.2 of |27|. One obtains the following - in leading order 
in Os and aem exact - expressions for the momentum-space wave functions of the virtual 



photon from (52) and (58) of 



,p'fl^{a,kT,Q) 



T —Qq 



± ((feO + m,)(A;'0 + nig) - a{l - a)q^) 6x,y6^^^. 



(A.l) 



+ e^^'l^'^kT \q\ 6x,^x'{a6^,±i - (1 - «)'^a,tO ± e^'^'^'k^^x,x'6x^^i 



1 



(4 + (feO + + a |q|)(/c'° + + (1 - a) \q\))5x,-x 
i(signA)</,fe^g-g^^)(^0 _ ^/o _ _ \q\)kT6x,X' 



+ e 



(A.2) 



Here q denotes the quark flavours, Qq their charges in units of the proton charge e = 
\/47raem- The number of colours is Nc = 3. The momenta k and k' of the quark and 



antiquark, respectively, are defined in (2^), and we have /ct = |^t| and (f^k = arg(fcTi + 
ikT2)- The normalisation factor is 



= e a{l - a) . 

Note that the wave functions satisfy the following relations 



Vig,^(a,fcT,Q) 



where 



i^^^l\a,kT,Q)] =^\\^,"(a,fc^,Q), 




(A.3) 

(A.4) 
(A.5) 

(A.6) 



With ( [A.4| ) and (|A.5| ) we can show that the distributions defined in ( p. 17 ) are real- 
valued. Indeed, from ( |2.12| ) we find that a^'^\kT, W) is independent of the direction of fc-r 
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and, therefore, 



a^'i\kT,w)j = a-('?)(-fcT,VF) 
= a-('?)(fc^,VF). 



(A.7) 



We then get from (12.17] ) 

d''aT,L{W,Q^V,v) 
drj df) 



Y.T, 

q A, A' 



d kr f d kx -j{q)±,L, r r/- a 



X fJ 



(27r)2 y (27r) 
(fcT - fcT, W)) * 6irj - AE) (M^'\a, kr, Q) 



q A, A' 



d^kx f d^kx ( 7(q)±,L 



(27r)2 J (27r)2 V "^AA- 
dW,L{W,Q^V,f]) 



X &(i\kf^ - k^, W) 5{ri - AE) V'iv^'^Ca, k^, Q) 



drj dr] 

In the high energy and small-ZcT hmit 



|q| ^ oo , Q < |q I , 
the photon wave functions simphfy to 



2 I 2 1 ~^ ^q 



< 1, 



k\^m\ 
(l-a)2q2 



« 1 



(A.8) 



(A.9) 



V;S,^(a,fcT,Q) 



^S^(a,fcT,Q) 



a(l - a)Q2 + ^2 + j„2 



"lg'^A,A"^A,±i 



± A:Te±''^'=5A,-A' (a^A.ii " (1 " «)'^A,Ti 



2^NcOLcraQq 



a(l — a)(52 _)_ /j2 _|_ TTij 



■ a(l - ajQ^A -A' • 



(A.IO) 
(A.ll) 



Inserting the expressions ( A.l(]| ), ( [A.llD into the Fourier transform ( |2.9| ) yields the standard 
formulae for the wave functions in longitudinal momentum and transverse position space 
(see again |27] for a detailed derivation) 



V^AA' 



(a,r,(5) = -^^yfa^Qq | 



A', -A 



'^A,±i - (1 - «) -^A.Ti ^l^li^l 



V^S^(a,r,Q) 



vr 



y/aeraQqQoi{l - a) Sy -\ Ko{eqr) . 



(A.12) 
(A.13) 



Here = arg(ri + ir2). 



eg = ^a(l - q)Q2 + ^2 ^ (A.14) 

and i^o,i are the modified Bessel functions. Hence we have in the high energy and small-fc^ 
limit ( [A.9| ) for the 'y*p cross sections the results ( |2.6D and (2.7), that is, 

^t,l(W^,Q') = V / dV /'daVUS,+'''(a,r-,Q) ' a^''Hr,W) (A.15) 
9 ^0 A,A' 
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Figure 5: Energy mismatch AE of the 7* dipole transition as a function of the absolute value 
/ct of the quark's transverse momentum (left graph) and of its longitudinal momentum fraction a 
(right graph). The parameters are W = 150 GeV, ^ iq QeV^, = 0.14 GeV. For the left 
graph we have chosen a = 1/2 and for the right graph kr — 1 GeV, respectively. The dotted line 
shows the absolute minimum Ai^niin.abs for the given external kinematics, that is, AE at kr — 
and a = 1/2. 



with 



Nr. 



!<,)+(«, r,Q) = ^ a^^Ql { [a' + (1 - af] el[K,{e,r)f + ml[Ko{e,r)f} , (A.16) 



A, A' 



^|V^if (a,r,Q) 

A,A' 



(A.17) 



For the calculation of ctt above we can choose either of the two polarisations, + or — in 
( A.15| ), since they lead to the same cross section. We have chosen the + polarisation here. 

Let us now discuss some basic properties of the energy mismatch AE ( 2.13| ) which 
follow directly from the kinematics of the 7* — > qq splitting. Similar considerations then 



apply to AE in (2.14). Explicitly we obtain from (2.13) and ( |2.8D for quark flavour q 
AE{kT, a) = Ja'^q^ + kl + ml + J {I - afq^ + + m? - g° . 



Figure |5| shows the dependence of AE{kT-,a) on the absolute value kx of the quark's 
transverse momentum (left graph) and on its longitudinal momentum fraction a for the case 
of light quarks (right graph). We note that AE{kT,a) strongly peaks at the longitudinal 
momentum endpoints and rises monotonically with the transverse momentum. At fixed a 
the energy mismatch becomes minimal for kx = 0: 



a^q^ + + /(I — af'q^ + w?„ — q . 



AE^M = AE{0,a) 



The absolute minimum of AE is reached at k^ = 0, a = 1/2, 



(A.19) 



A^niin,abs = AE(0, 1/2) = J + Am^ - . 



(A.20) 
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That is, we have 

AE{kT,a) > AEmin(a) > A£;^in,abs > , (A.21) 

where the last inequality is strict for all > 0. Thus we see that there is an a priori 
minimal value for the energy mismatch AE at a given kinematical point and for given 
quark flavour. Figure ^ in section |3| shows twice this minimal energy mismatch. 

B. Some technical details of the calculation 

In this appendix we give the details of the calculations for the AE distributions ( 2.17] ), 



( 2.19 ), and the corresponding quantities for individual quark flavours q. We use the Golec- 



Biernat-Wiisthoff model pT| which describes the F2 data from HERA quite well. The 



dipole cross section of this model is given in (3.1) and ( |3.2[ ). We point out again that 
the dipole cross sections depend not only on W but also on through the x 

dependence. We use the GBW model nevertheless, since we shall only be concerned with 
specific kinematic values of W and for which we study loffe-time distributions. We 
shall not compare structure functions at the same W and different values, where the 



choice of energy variable is essential, see |28, p9[. 



The Fourier transformation ( 2.12 ) of d^g^(r, x) ( |3.1| ) gives 

~a^^ly^{kT, x) = (70 [(27r)2<5(2) (kr) - Airr^o exp(-r24)] . (B.l) 

Our aim is to calculate the AE^ distributions for the various quark-flavour contributions 
to ax and aL- We have 



(B.2) 



X (^W(fcT - kr, W) 6iv - AE) M^'\a, kr, Q) 



da^^\{W,Q\AE, 
dAE. 











dri drj 



We insert the photon wave functions from ( A.l ) and ( A. 2 ) and replace a^'^\kT — kx, W) by 
^GBW^^T — kx, x) from (|B.1| ). In order to integrate out the azimuthal angles we decompose 
the photon wave functions as follows 



V;S^'^(a,fcT,Q)= ^ e-<^'=^g,^'^(a,feT,Q). 



(B.4) 



n=-2 



For any function /(/cr) we have the relation 



Jo 



f{kT{AE,a)). 



(B.5) 
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Using this we find for the distributions in AE^ that 



da. 



(g) 

T,L 



da. 



(g), const 
T,L 



da. 



dAE. 



dAE. 



+ 



(g),e 

T,L 



dAE. 



with 



da. 



(g), const 



dAE4 



iAE+) 

■ao I da 



da. 



T,L 



dAE. 



-ao 



{AE+) fAEn,i^{a) 

da / dAE. 



\,\'n 



(B.6) 



(B.7) 



t,0L/0i,0i,/0„2 

r\j nj hj r\j / q 



In{2rlkTkT) J{q) 



-rl{kT-~kT? 



A,A',n 



f exp(2r2/cyA;7.) 



(B. 



where are the modified Bessel functions. Here, the energies and momenta of quark and 
anti-quark (^.§|) are understood as functions of a, AE, AE via 



kT 



kT 



'{{AE + g0)2 - g2) ((A£; + g0)2 _ (1 _ 2a)2q2) 



4(A£; + 



o^2 



((A^ + q0)2 - q2) ((A^ + q0)2 _ (1 _ 2a)2q2) 



4(AE + gO)2 



m„ 



(B.9) 



(B.IO) 



The a-integration range in ( p.T] ) is finite for fixed A£^+ and its endpoints correspond to 
vanishing transverse momenta. The extremal values of a are given by the two solutions of 
the equation AS_|_/2 = A£^|^^^g with respect to a: 



a 



max, mm 



2 2|g| V 



4m2 



{AE+/2 + q' 



0\2 



i2 ■ 



(B.ll) 



In all cases considered in this paper amin ~ and amax ~ 1- 

We perform the residual integrations in ( [B.7D and (|B.8| ) numerically. We avoid nu- 
merical errors due to the large cancellations in sums like ( |A.1§| ) by proper rewriting, and 
we make sure that the integration error is under control by using different integration 



algorithms pq, |36|. 
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